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Abstract 


Measurements  of  the  luminescence  lifetime  of  Ru(bpy)3^* 

(bpy  =  2 , 2 ' -bipyridine)  in  Nafion  were  used  to  probe  the 
heterogeneity  of  -SO3H  sites  in  the  Nafion  membrane.  The  decay 
kinetics  of  the  probe  ion  are  shown  to  depart  from  simple  first- 
order  behavior.  Albery's  model  of  a  continuous  distribution  of 
rate  constants  is  used  to  explain  the  decay  kinetics.  The 
results  of  these  and  related  investigations  indicate  that  the 
probe  ions  reside  in  sites  with  different  local  water  activities 
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Introduction 


Ion-containing  polymers  (ionomers)  have  a  wide  variety  of 
applications  including  uses  in  electrochemistry,  chemical 
synthesis,  and  separation  science. Nafion,  a  family  of 
perfluorosulfonate  ionomers  marketed  by  Du  Pont,  has  proven  to  be 
an  extremely  useful  and  versatile  ion-containing  polymer.  Nafion 
and  other  ionomers  have  an  unusual  morphology  in  which  the 
charged  groups  aggregate  to  form  domains  called  ionic  clusters. 
These  clusters  resemble  reversed  micelles  and  are  randomly 
distributed  throughout  the  backbone  chain  material  phase.  The 

precise  details  of  this  structure  are  still  unknown. 

Information  on  the  structural  and  dynamic  nature  of  Nafion 
is  a  prerequisite  for  the  future  utilization  of  this  material. 
Specifically,  information  about  the  chemical  microenvironment  and 
diffusion  of  ions  and  molecules  within  the  ionic  clusters  is 
required.  Luminescence  probe  techniques  can  provide 
photophysical  and  photochemical  information  about  the  properties 
of  this  system.  For  this  reason,  we  and  others^®'^^  have  been 
using  luminescence  probe  techniques  to  study  the  morphology  and 
chemical  characteristics  of  the  ionic  cluster  phases  in  ionomers. 

The  specific  luminescent  probe  molecule  used  in  this  study 
is  tris(2,2'-bipyridyl)  ruthenium ( II )  ,  Ru(bpy)3^'^,  an  extensivf-ly 
studied  metal  complex.^®'®®  Lee  and  Meisel  reported  luminescence 
quenching  studies  of  Ru(bpy)3^^  in  Nafion.^®  Lee  and  Meisel 
observed  that  the  luminescence  decay  of  Ru(bpy)3^^  in  Nafion  was 
monoexponential;  in  contrast,  we  have  always  obtained 
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nonexponential  decays  for  the  luminescence  of  Ru(bpy)3^*  in 
Nafion.  Possible  reasons  for  the  discrepancies  between  our 
results  and  those  of  Lee  and  Meisel  will  be  presented  here. 

The  nonexponential  luminescence  decay  curves  reported  in 
this  paper  can  be  fitted  using  a  biexponential  decay  model. 
However,  by  fitting  the  data  to  this  model  it  is  tacitly  assumed 
that  Ru(bpy)3^*  exists  in  two  distinct  chemical  environments 
within  the  Nafion  membrane.  An  alternative  to  fitting  the  data 
to  the  biexponential  model  is  to  assume  that  a  continuous 
distribution  of  chemical  microenvironments  exists  within  the 
membrane.  This  distribution  of  microenvironments  produces  a 
continuous  distribution  of  luminescence  decay  rate  constants. 

We  have  used  Albery's  model  for  dispersed  kinetics  to 
analyze  the  luminescence  of  Ru(bpy)3^'^  in  Nafion  membrane. The 
fits  of  the  experimental  data  to  the  dispersed  kinetics  model  are 
better  than  the  fits  to  the  biexponential  decay  model.  This  data 
analysis  shows  that  the  decay  of  Ru(bpy)3^*  in  Nafion  can  be 
described  by  a  Gaussian  distribution  of  the  energy  of  activation 
for  the  decay.  Information  about  the  chemical  distribution  and 
nature  of  the  ion  exchange  sites  in  the  polymer  phase  can  be 
obtained  from  this  analysis.  These  and  related  results  are 
presented  in  this  paper. 

Experimental  Section 

MATERIALS .  Nafion  117  membrane  (proton  form)  was  obtained  from 
the  E.I.  Du  Pont  Company.  Ru(bpy)3Cl2  SHjO  was  obtained  from  G. 

F.  Smith  and  was  used  as  received.  Water  was  circulated  through 
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a  Milli-Q  water  purification  system  (Millipore  Corp.)  or  was 
pyrolytically  triply  distilled.  All  other  reagents  and  solvents 
were  of  the  highest  available  grade  and  ware  used  without  further 
purification.  All  glassware  was  cleaned  in  2M  HNO3  which  had 
been  doubled  distilled  in  quartz  (Seastart) . 

PROCEDURES .  As-received  Nafion  membrane  was  cut  into  ca.  10  X  30 
mm  rectangles.  The  precut  membranes  were  then  cleaned  by 
ultrasonicating  in  several  portions  of  pyrolytically  triply 
distilled  water  (PTD-H2O)  ,  and  then  ultrasonicating  in  50:50 
EtOH-PTD-HjO  (0.5  hr.).  The  membranes  were  then  boiled  in  PTD- 
HjO  for  2  hr.  and  rinsed  ultrasonically  in  PTD-H2O.  The 
membranes  were  changed  to  the  Na"^  form  by  soaking  overnight  in  IM 
NaOH,  and  finally  rinsed  ultrasonically  in  PTD-H2O. 

The  wet  membranes  were  weighed  separately  and  placed  in 
glass  vials.  Two  of  the  membranes  were  used  to  determine  the 
equilibrium  membrane  water  content.  These  two  membranes  were 
dried  in  a  vacuum  oven  overnight  at  115°C  and  reweighed.  The 
water  content  was  calculated  as  percent  by  weight  (100  X  (weight 
of  H2O)  /  (weight  of  dry  polymer) ) . 

The  luminescent  probe  Ru(bpy)3^'^  was  incorporated  (loaded) 
into  the  remaining  Nafion  membranes  by  stirring  each  membrane  in 
a  separate  vial  containing  the  required  concentration  of 
Ru(bpy)3^*.  (The  aqueous  solution  of  Ru(bpy)3^*  used  to  load 
Ru(bpy)3^*  into  the  membranes  is  called  the  "loading  solution"; 
each  membranes  was  placed  in  a  loading  solution  containing  the 
precise  concentration  of  Ru(bpy)3^*  needed  to  load  the  membrane  to 
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the  desired  Ru(bpy)3^^  level.)  As  will  be  discussed  in  detail 
later,  because  diffusion  of  the  probe  into  the  polymer  phase  is 
slow,  the  membranes  were  equilibrated  with  their  respective 
loading  solution  for  1  month  before  use.  The  quantity  of 
Ru(bpy)3^*  loaded  was  maintained  very  low;  between  0.8  and  0.01% 
of  the  ionic  sites  in  the  Nafion  membranes  were  exchanged  with 
the  luminescent  probe.  After  loading,  and  previous  to 
measurement,  the  membranes  were  rinsed  with  water. 

Some  experiments  were  performed  on  oxygen-free  Nafion 
membranes.  Oxygen  was  removed  from  the  Nafion  membranes  as 
follows:  After  loading  with  Ru(bpy)3^^  the  membrane  was  placed  in 

a  12.5  mm  X  12.5  mm  X  45  mm  rectangular  cuvette  filled  with 
water.  The  cuvette  was  sealed  with  a  rubber  septum,  and  nitrogen 
was  bubbled  into  the  solution  for  one  and  a  half  hours  while  the 
solution  was  ultrasonicated. 

SPECTROSCOPY .  Steady-state  emission  spectra  were  obtained  with  a 
Spex  Fluorolog  2  spectrof luorometer .  The  samples  were  excited 
with  a  450-W  xenon  lamp  and  all  the  slits  were  set  at  1.25  mm. 

The  excitation  wavelenght  was  452  nm.  Luminescence  was  detected 
perpendicular  to  the  incident  radiation.  Luminescence  lifetime 
measurements  were  obtained  using  a  time-correlated  single  photon 
counting  spectrometer  (Edinburgh  Instruments  Model  199) .  The 
spectrometer  is  fitted  with  a  nitrogen  flashlamp  (337  nm)  and  a 
microcomputer  (IBM  PS/2  Model  30  with  a  8087  microprocessor) . 

The  Ru(bpy)  3^^- loaded  membranes  were  sandwiched  between  two 
quartz  slides.  This  assembly  was  positioned  in  the  spectrometer 
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at  45  degrees  to  the  axis  of  excitation.  Emission  was  detected 
from  the  opposite  side  of  the  membrane  at  45  degrees  to  the  axis 
of  the  detector.  This  orientation  (and  the  emission  filter) 
minimizes  the  detection  of  scattering  from  the  membranes. 

Data  Analyses.  The  general  approach  for  obtaining  kinetic 
parameters  from  the  luminescence  decay  curves  was  described 
previously. We  fitted  the  experimental  intensity  vs.  time 
transients  to  specific  decay  models  using  the  method  of  nonlinear 
least-squares^^'^*  (Marquardt  algorithm)  ;  the  Edinburgh 
Instrument's  or  MTR  Software  Inc.'s  nonlinear  least-squares 
software  was  used.  The  Marquardt  algorithm  minimizes  the 
goodness-of-f it  parameter  chi  square  (x^)  given  by 

=  S  {  [y^  -  y(ti)  f)  (1) 

where  are  the  uncertainties  (errors)  in  the  data  points  yi  and 
y(ti)  is  the  value  of  the  fitting  function  at  each  data  point. 

For  single  photon  counting  data 

Three  different  decay  models  were  used  to  fit  the  decay 
curves.  The  first  two  models  were  either  a  monoexponential  or  a 
biexponential  decay  model. These  models  assume  that  one  or  two 
populations  of  emitters  exist  in  the  system. 

The  third  model  used  to  fit  our  data  is  Albery's  dispersed 
kinetics  model,  which  takes  into  account  the  heterogeneous  nature 
of  the  system.  The  luminescence  of  molecules  emitting  from 
heterogeneous  systems  (e.g,  proteins,  surfaces,  and  micelles) 
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Therefore,  a  plot 


frequently  d='part  from  first-order  kinetics. 
of  the  loga^^thm  of  emission  intensity  vs.  time  is  nonlinear. 
Albery's  model  accounts  for  this  departion  from  first-order 
kinetics. 

Albery's  model  is  an  adaptation  of  a  new  model  for 
interpreting  the  decay  kinetics  in  nonhomogeneous  systems . 

Rather  than  assuming  just  one  or  two  distinct  populations  of 
emitters,  this  model  assumes  that  a  continuous  distribution  of 
decay  rates  exits  within  the  heterogeneous  system.  This 
continuous  distribution  of  decay  rates  results  from  a  continuous 
distribution  of  chemical  microenvironments  within  the  system; 
each  microenvironment  has  a  distinct  decay  rate  constant.  This 
model  assumes  that  each  of  these  distinct  rate  constants  is 
first-order. 

Recently,  Albery  et  al.®^  proposed  this  general  model  for 
dispersed  kinetics  in  heterogeneous  systems.  This  model  assumes 
that  the  dispersion  in  the  kinetics  is  the  result  of  a  change  in 
energy  of  activation  for  the  reaction  at  different  sites  in  the 
system.  Albery's  model  assumes  that  the  observed  kinetics  can  be 
described  by  a  Gaussian  distribution  of  the  energies  of 
activation,  or  of  the  natural  logarithm  of  the  decay  rate 
constant. 

We  have  used  Albery's  model  as  the  third  approach  for 
fitting  our  experimental  emission  decay  curves.  The  equation 
used  to  fit  the  decay  curves  to  this  Gaussian  distribution  model 


C/Co  = 


I  exp(-x^)  exp[-r  exp(7x)  ]  dx 

p-HX) 

J  exp(-x^)  dx 


(2) 


where 


i: 


exp(-x  )  dx  =  TT 


-  _l/2 


C  is  the  concentration  of  excited  state  species  at  any  time,  t, 
after  excitation,  is  the  excited  state  concentration  at  t  =  0, 
X  is  the  parameter  that  describes  the  Gaussian  (see  Appendix) ,  7 
corresponds  to  the  width  of  the  distribution,  and  r  =  )ct.  The 
decay  rate  constant  of  maximum  probability  is  Ic.  The  r  in 
Equation  4  should  not  be  confused  with  the  average  lifetime  for 
the  emitter  which  is  1/Ic.  The  integral  in  Equation  2  is  solved 
numerically  as  described  by  Albery  et  al.®^  Equation  2  is 
derived  in  the  Appendix. 

Albery 's  dispersed  kinetics  model  provides  the  "most 
probable”  rate  constant,  k,  which  is  the  first-order  rate 
constant  associated  with  the  average  chemical  microenvironment 
within  the  medium.  This  rate  constant  can  be  studied  in  the 
normal  fashion  as  if  the  process  was  exponential  (i.e.,  normal 
kinetic  analysis  can  be  performed  with  this  "most  probable"  rate 
constant  value) .  The  dispersion  of  k  (the  7  value)  provides  a 
measure  of  the  disorder  of  the  medium.  This  information  can  be 
used  to  investigate  the  local  structure  of  the  medium  and  the 
dependence  of  the  rate  process  on  energetic  and  structural 
parameters.*® 
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Results  and  Discussion 

Previous  Investigation  of  RufbPvK^^  Emission  from  Nafion .  Lee 
and  Meisel^®  observed  monoexponential  decay  curves  for  the 
emission  of  Ru(bpy)3'‘  in  Nafion.  We  have  repeated  this 
experiment  numerous  times  over  the  last  5  years  and  have  never 
observed  monoexponential  decays.  Recently,  Kaneko  and  Hayakawa'’^ 
also  observed  that  the  decay  of  Ru(bpy)3^*  in  Nafion  films  cannot 
be  describe  by  a  monoexponential  decay  model.  Figure  1  shows 
typical  decay  curves  for  Ru(bpy)3^‘^  in  Nafion  at  various  loading 
levels;  note  that  these  are  semi  logarithmic  plots.  The 
nonlinearity  observed  i.n  Figure  1  clearly  shows  that  a 
monoexponential  decay  model  is  not  appropriate.^^  Nonlinear 
semilogarithmic  plots  such  as  these  have  been  observed 
consistently  by  our  group,  independent  of  instrument,  operator, 
or  time. 

We  believe  that  it  is  possible  that  Lee  and  Meisel's 
results  were  affected  by  a  nonequilibrium  distribution  of 
luminescent  probe  ions  in  the  Nafion  membrane.  Lee  and  Meisel 
stirred  their  membranes  in  the  Ru(bpy)3^'^  loading  solution  for 
only  12  hours.  A  simple  calculation  shows  that  equilibrium  times 
well  in  excess  of  12  hours  are  required  if  probe-membrane 
equilibrium  is  to  be  achieved.  The  distance,  a,  a  diffusing 
species  moves  in  time,  t,  is  given  by  A  =  (2Dt)^^^  where  D  is  the 
diffusion  coefficient.^®  Assuming  D  for  Ru(bpy)3^'^  in  Nafion  is 
2  X  10'^°  cmVs,^  a  Ru(bpy)3^‘^  ion  moves,  on  average,  only  0.004  cm 
into  a  Nafion  membrane  during  a  12  hour  equilibrium  period.  In 
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order  to  achieve  equilibrium  the  ion  must  move  half-way  through 
the  membrane  (0.0125  cm).  Thus,  after  12  hours,  a  nonhomogeneous 
distribution  of  Ru(bpy)3^*  is  extant  in  the  membrane;  electron 
microprobe  analyses  have  experimentally  confirmed  such 
nonhomogeneous  distributions.^^ 

It  is  not  obvious  why  an  nonhomogeneous  distribution  of 
Ru(bpy)3^'^  in  the  membrane  would  yield  a  monoexponential  decay. 

We  have  found,  however,  that  after  loading,  the  wavelength  of 
maximum  emission  intensity  for  Ru(bpy)3^‘^  blueshifts  for  periods 
of  up  to  one  month.  This  blue  shift  results  from  a  gradual 
partitioning  of  the  probe  from  domains  with  high  local  water 
activity  to  domains  with  high  local  fluorocarbon  chain  material 
activity  (vide  infra  and  ref.  19) .  These  data  suggest  that 
Ru(bpy)3^'^,  in  freshly-loaded  films,  is  predominately  present  in 
an  aqueous-like  environment.  Since  only  this  aqueous-like 
environment  is  populated  at  short  equilibration  times,  a 
monoexponential  decay  (as  observed  by  Lee  and  Meisel)  might  be 
anticipated.  In  contrast,  for  long  equilibration  periods  a 
variety  of  microenvironments  become  populated  (aqueous-like  and 
fluorocarbon-like)  and  this  produces  the  nonexponential  decay 
observed  in  our  experiments. 

Biexponential  Decay  Model.  The  biexponential  decay  model  assumes 
that  there  are  two  distinct  populations  of  emitters  (or  two 
distinct  chemical  microenvironments)  within  the  membrane  phase. 

As  a  result  of  these  populations,  two  distinct  decay  rate 
constant  or  decay  lifetimes  are  observed.  Figure  1  shows  that 
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our  data  can  be  fitted  to  a  biexponential  decay  model.  The 
lifetimes,  standard  deviations,  and  values  obtained  from  these 
fits  are  shown  in  Table  I. 

The  kinetic  data  obtained  from  the  biexponential  decay  model 
(Table  I)  may  be  summarized  as  follows.  At  loading  levels  less 
than  0.07%,  both  the  long  and  the  fast  lifetimes  increase  with 
the  quantity  of  Ru(bpy)3^*  loaded  into  the  membrane.  Furthermore, 
the  data  in  Table  I  suggest  that  the  microenvironment  yielding 
the  long  lived  species  is  preferentially  occupied.  (At  low 
loading  levels  84  to  98%  of  the  emitters  occupy  sites  in  this 
microenvironment,  while  at  higher  loading  levels,  only  68  to  84% 
of  the  emitters  are  in  this  long  lifetime  environment.)  Above 
0.07%  loading,  the  lifetimes  and  relative  sizes  of  the  two 
populations  remain  more  or  less  constant  (Table  I) .  Kaneko  and 
Hayakawa*®  described  the  decay  of  Ru(bpy)3^'^  adsorbed  in  Nafion  at 
a  low  loding  level  as  a  biexponential  decay  with  long  lifetime 
component  of  750-839  ns  (over  90%)  and  a  fast  lifetime  component 
of  121-238  ns. 

The  kinetic  data  in  Table  I  suggest  that  the  biexponential 
decay  model  is  not  appropriate  for  this  system.  The  first 
problem  with  this  model  concerns  the  changes  in  both  the  short 
and  long  lifetimes,  observed  at  low  loading  levels.  If  two 
distinct  populations  exist  in  the  membrane,  the  two  lifetimes 
obtained  from  the  kinetic  analysis  should  be  the  same  regardless 
of  the  loading  level;  the  relative  numbers  of  emitters  in  each 
population  might  change  with  loading,  but  there  is  no  a  priori 
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reason  to  expect  that  the  lifetime  will.  The  change  in  lifetimes 
observed  in  Table  I  suggest  that  more  than  two  distinct 
environments  are  being  populated. 

The  data  obtained  from  the  biexponential  model  also 
contradicts  evidence  obtained  from  steady-state  emission 
experiments.  Table  II  shows  that  the  wavelength  of  maximum 
emission  intensity  (A„„)  is  strongly  blue  shifted  (relative  to 
emission  from  water)  at  low  loading  levels;  shifts  to  more 

aqueous-like  values  at  higher  loading  levels.  The  blue-shifted 
clearly  shows  that  the  first  Ru(bpy)3^^  ions  which  enter  the 
membrane  occupy  sites  with  high  local  fluorocarbon  chain  material 
activity  (and  low  local  water  activity) When  these  sites  are 
occupied  the  subsequently-loaded  Ru(bpy)3^*'s  enter  sites  with 
higher  local  water  activities;  this  produces  the  shift  in  A„„  to 
more  aqueous-like  (and  less  fluorocarbon-like“)  values. 

While  the  lifetime  data  obtained  from  the  biexponential 
decay  model  (Table  I)  show  evidence  for  preferential  occupancy  of 
one  environment  over  the  other,  the  lifetimes  are  reversed.  Note 
that  the  results  from  the  biexponential  model  (Table  I)  suggest 
that  at  low  loading  the  long  (aqueous-like)  lifetime  is 
preferentially  occupied  (Table  I) .  This  is  in  direct 
contradiction  to  the  steady-state  data  which  suggest  that  the 
nonaqueous  or  fluorocarbon-like  environment  is  preferentially 
occupied  (Table  II) . 

The  final  problem  with  the  biexponential  decay  model  data 
concerns  the  very  short  lifetime  observed  at  low  loading  levels. 
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It  is  difficult  to  explain  how  such  a  low  lifetime  for  Ru(bpy)3^"^ 
could  be  physically  achieved;  for  example,  while  Ru(bpy)3^^'s 
emission  parameters  from  numerous  solvents  have  been 
investigated,^®  none  of  these  solvents  yield  such  short 
lifetimes.  One  possible  explanation,  which  will  be  explored  in 
detail  below,  is  that  the  lifetime,  at  low  loading,  is  decreased 
by  O2  quenching.  However,  the  concentration  of  O2  in  the 
membrane^®  is  too  low  to  yield  such  dramatically  reduced 
lifetimes. 

The  bottom  line  is  that,  while  it  is  impossible  to  state 
with  100%  assurance  that  a  model  is,  or  is  not,  appropriate  for  a 
particular  system,  the  bulk  of  the  evidence  suggests  that  the 
biexponential  decay  model  is  not  appropriate  here.  Furthermore, 
as  will  be  shown  below,  the  dispersed  kinetics  model  provides 
statistically  better  fits  to  the  experimental  data. 

Dispersed  Kinetics  Model.  The  data  in  Tables  I  and  II  suggest 
that  more  than  two  chemical  microenvironments  are  present  in 
Nafion  membrane;  this  has  been  confirmed  by  Mossbauer 
investigations  of  Fe®*-exchanged  Nafion®*’®®  and  infrared  and  NMR 
studies  of  water  in  Nafion  membranes.®®’®^  Faulkner  and  Bartolo 
have  recently  proposed  that  a  continuous  distribution  of  chemical 
microenvironments  exists,  for  Ru(bpy)3®^,  within  the  cation 
exchange  polymer  poly (styrenesulfonate)  (PSS).®®  These  authors 
used  the  dispersed  kinetics  model  of  Albery  et  al.®^  to  analyze 
Ru(bpy)3®^  emission  data  for  PSS. 
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Figure  2  shows  that  Albery's  dispersed  kinetics  model  can  be 
fitted  to  the  experimental  Ru(bpy)3^^  decay  data  obtained  here. 

Two  independently  adjustable  parameters  k,  the  average  decay  rate 
constant  and  7,  the  width  of  the  distribution,  were  used  to 
obtain  these  fits.  These  kinetic  parameters,  as  well  as  the 
values  for  the  fits,  are  shown  in  Table  III. 

The  values  obtained  using  the  dispersed  kinetics  model 
(Table  III)  are  essentially  identical,  at  all  loading  levels,  to 
the  values  obtained  from  the  biexponential  model  (Table  I) . 
This  would  seem  to  suggest  that  neither  model  is  preferable; 
however,  four  independently  adjustable  parameters  (Two  k's  and 
two  B's)  were  used  in  the  biexponential  fit,  whereas  only  two 
adjustable  parameters  were  used  for  the  dispersed  kinetics  model. 
Thus,  in  addition  to  the  arguments  presented  above,  the  principle 
of  "Ockham's  Razor"^®  makes  the  dispersed  kinetics  model 
preferable  to  the  biexponential  model.  More  important  than  the 
philosophical  "Ockham's  Razor"  argument,  the  standard  deviation 
for  the  parameters  obtained  from  the  fit  to  the  dispersed 
kinetics  model  were  always  smaller  than  for  the  biexponential  fit 
(Table  I  and  III) .  This  provides  guantitative  statistical 
evidence  that  the  dispersed  kinetics  model  provides  a  better  fit 
to  the  experimental  data. 

Albery's  dispersed  kinetics  model  is  a  functional 
mathematical  model  that  can  be  used  to  obtain  kinetic 
information.  How  can  this  information  be  translated  into  a 
chemical  model?  A  chemical  picture  can  be  envisioned  in  which 
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the  clusters  of  Nation  offer  a  variety  of  chemical 
microenvironments  to  the  Ru(bpy)3^^  ions.  These  microenvironments 
arise  because  different  parts  of  the  clusters  have  different 
dielectric  constants,  depending  on  the  amount  of  water  and 
per fluorocarbon  chain  material  present. 

The  characteristics  of  a  specific  chemical  microenvironment 
can  have  a  significant  impact  on  the  Ru(bpy)3^*  decay  rate.  The 
lifetiem  of  Ru(bpy)3^^'s  metal-to-ligand  charge  transfer  (MLCT) 
excited-state  decay  is  dependent  on  both  radiative  and 
nonradiative  modes  of  deactivation. In  addition,  the  excited- 
state  decay  is  dependent  on  thermal  activation  to,  and  subsequent 
decay  from,  low-lying  dd  states.  These  processes  are  known  to  be 
dependent  on  the  chemical  microenvironment  in  which  Ru(bpy)3^'^ 
resides. Thus,  the  excited-state  lifetime  will  vary  from 
microenvironment  to  microenvironment  within  the  membrane. 

Recent  studies  have  demonstrated  the  effects  of  the 
properties  of  the  chemical  microenvironment  on  the  emission 
characteristics  of  Ru(bpy)3^^.  Meyer  et  al.®°  have  observed 
nonexponential  decays  for  a  Ru(II)  polypyridyl  complex  bound  to  a 
chlorosulfonated  polystyrene  film.  The  observed  lifetimes  were 
dependent  on  the  hydrophobic ity  and  structure  of  the  chemical 
sites  within  the  film.  Furthermore,  Vining  and  Meyer®^  have 
presented  electrochemical  and  spectroelectrochemical  evidence  for 
at  least  three  different  chemical  regions  within  Nafion  films. 

The  presence  of  oxygen  in  these  regions  brings  another  decay 
mechanism  for  the  deactivation  of  *Ru(bpy)3^*. 


14 


The  k  values  obtained  from  the  fits  to  the  dispersed  kinetic 
model  can  be  used  to  calculate  the  mean  lifetime,  r,  values  for 
each  Ru(bpy)3^*  loading  level.  These  r  values  are  shown  in  Table 
III.  The  r  value  for  the  lowest  loading  level  appears  to  be 
significantly  lower  than  the  other  r's.  Applying  the  standard  Q 
test  for  rejection  of  data,®^  we  find  that  at  the  99%  confidence 
level,  the  r  associated  with  the  lowest  loading  level  is 
significantly  less  than  the  r's  obtained  at  the  high  loading 
levels.  Furthermore,  the  average  r  value  obtained  from  the 
higher  (0.03  to  0.80%)  loading  levels  (435  ±  8  ns)  is  identical 
to  the  lifetime  for  Ru(bpy)3^*  in  aerated  aqueous  solution  (434 
ns)  . 

The  trends  in  the  kinetic  parameters  obtained  from  the 
dispersed  kinetics  model  (Table  III)  agree,  phenomenologically, 
with  the  trends  in  the  steady-state  emission  data  (Table  II) . 

Both  data  show  that  an  aqueous-like  environment  is  populated  at 
high  loading  levels  whereas  a  "non-aqueous"  environment  is 
preferentially  populated  at  low  loading  levels.  As  noted  above, 
this  "non-aqueous"  environment  simply  means  that  the  first 
Ru(bpy)3^*'s  which  enter  the  membrane  preferentially  occupy  sites 
with  high  local  fluorocarbon  chain  material  activities  and  low 
water  activities;  given  the  known  hydrophobicity  of  Ru(bpy)3^*, 
this  preference  for  hydrophobic  sites  is  not  surprising.®^ 

The  width  of  the  distribution  (y  values  in  Table  III)  also 
follow  the  trends  observed  in  the  r  and  values;  as  the 
loading  level  is  increased  the  overall  Gaussian  distribution  is 


15 


reduced  in  width.  This  reduction  in  width  occurs  because  many  of 
the  component  lifetimes  are  similar  in  the  aqueous-like  domains 
and  because  a  relatively  smaller  number  of  the  fast-decaying 
species  are  present.  Therefore,  at  higher  loading  levels  there 
is  a  smaller  contribution  of  the  ions  in  the  "non-aqueous”  sites 
to  the  total  distribution. 

The  question  which  arises,  however,  is  why  would  a  high 
local  chain  material-activity  microenvironment  yield  a  diminished 
lifetime  for  Ru(bpy)3^*.  We  believe  that  part  of  this  diminution 
in  Ru(bpy)3^^  lifetime  results  from  a  higher  local  O2  activity. 

It  is  well  known  that  Oj  is  more  soluble  in  fluorocarbon- 
containing  media  than  in  pure  water. For  example,  Lee  and 
Rodgers®^  have  shown  that  singlet  oxygen  is  five  times  more 
soluble  in  the  perfluorocarbon  region  than  in  the  water  clusters. 
Thus,  we  propose  that  the  "non-aqueous”  sites  are  also  sites  of 
high  local  Oj  activity.  Thus,  the  Ru(bpy)3^*'s  which  occupy  these 
high  fluorocarbon  chain  material-activity  sites  have  shorter 
lifetimes  due  to  the  higher  local  quencher  concentrations. 

The  effect  of  the  local  Oj  activity  on  the  lifetime  of 
Ru(bpy)3^*  was  confirmed  by  performing  luminescence  lifetimes 
studies  of  Ru(bpy)3^*  in  Nafion  membranes  that  have  been  deaerated 
(oxygen  removed) .  Figure  3  compares  two  luminescence  decay 
curves  for  the  same  loading  level  of  Ru(bpy)3^^  in  aerated  and 
deaerated  Nafion  membranes.  Figure  3  clearly  shows  that  the 
luminescence  decay  for  the  deaerated  membrane  is  longer  lived 
than  the  luminescence  decay  in  the  aerated  membrane. 
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Table  IV  shows  the  kinetic  parameters  obtained  from  the  fit 
of  the  luminescence  decay  curves  from  the  deaerated  membranes  to 
the  dispersed  kinetics  model.  In  agreement  with  the  data  in 
Figure  3,  the  lifetimes  for  the  deaerated  membranes  are  longer 
than  the  corresponding  lifetimes  in  the  aerated  membranes.  This 
observation  confirms  that  oxygen  quenches  the  luminescence  of 
Ru(bpy)3^*  in  the  aerated  membranes. 

The  kinetic  parameters  for  the  deaerated  membranes  (Table 
IV)  show  the  same  trends  as  the  kinetic  parameters  for  the 
aerated  membranes  (Table  III).  The  Ru(bpy)3^*  luminescence 
lifetime  increases  at  higher  loading  levels  for  both  the 
deaerated  and  aerated  membranes.  The  increase  in  luminescence 
lifetime  indicates  that  the  later-loaded  Ru(bpy)3^*  ions  occupy 
different  chemical  microenvironments  than  the  first-loaded  metal 
ions.  These  results  prove  that  a  multiplicity  of  chemically 
distinct  microenvironments  exists  in  Nafion  irrespectively  of  the 
presence  of  Oj.  In  addition.  Table  IV  shows  that  the  width  of 
the  distribution  (7  value)  is  reduced  as  the  Ru(bpy)3^^  loading 
increases  in  the  deaerated  membranes.  As  discussed  for  the 
aerated  membranes,  this  reduction  in  the  width  of  the 
distribution  occurs  because  the  relative  number  of  fast  decaying 
species  is  reduced  at  these  high  loading  levels. 

Previous  studies^®’^®  have  suggested  that  Ru(bpy)3^^ 
preferentially  resides  at  the  interface  between  the  water  cluster 
and  the  per fluorocarbon  backbone  region.  Furthermore,  Lee  and 
Meisel^^  observed  that  oxygen  quenches  the  fluorescence  of  pyrene 


17 


in  Nafion;  pyrene  also  resides  at,  or  near,  the 
fluorocarbon/water  cluster  interphase. The  results 
presented  here,  and  Lee  and  Meisel's  results, indicate  that 
oxygen  can  access  the  water  clusters  in  Nafion  to  quench  the 
pyrene  or  the  Ru(bpy)3^^  luminescence. 

Finally,  Table  IV  shows  that,  at  high  loading  levels,  the 
width  of  the  distribution  is  smaller  in  deaerated  membranes  than 
in  aerated  membranes  (Table  III) .  The  smaller  width  of  the 
distribution  occurs  because  the  absence  of  O2  reduces  the  number 
of  chemically  distinct  sites  in  Nafion.  This  is  in  agreement 
with  Lakowicz  et  al.''°  who  proved  that  quenching  causes  a 
lifetime  distribution  to  become  wider  and  shift  toward  shorter 
decay  times. 

Conclusions 

The  dispersed  kinetics  model  used  here  assumes  a  Gaussian 
distribution  of  the  logarithm  of  the  decay  rate  constant. 

However,  other  distribution  models  are  possible  (e.g., 

Lorentzian,  which  generally  gives  narrower  distributions).*®  The 
distribution  model  used  in  the  fitting  equations  could  have  also 
been  a  Gaussian  or  Lorentzian  distribution  of  rate  constants,  or 
lifetimes,  or  energies  of  activation,  etc.  In  addition,  there  is 
no  reason  to  believe  that  the  "real"  lifetime  distribution  is 
symmetric.  Finally,  it  is  possible  that  a  second  distribution 
component  should  be  included,  making  the  distribution  bimodal. 

At  present  there  is  no  theoretical  model  to  predict  the  exact 
distribution  and  therefore  the  selection  of  a  particular 
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distribution  model  is  made  arbitrarily.  However,  the  key 
question  is  whether  the  decay  process  is  best  described  by 
discrete  lifetimes  or  a  continuous  distribution  of  lifetimes;  in 
the  present  case  a  distribution  model  is  appropiate. 

The  results  obtained  here  prove  that  oxygen  is  not  the  sole 
reason  for  the  heterogeneity  in  the  luminescence  decay  kinetics 
observed  in  Nafion.  These  results  provide  evidence  for  the 
existence  of  chemically  distinct  regions  within  Nafion  which 
affect  the  lifetime  of  the  excited  state  of  Ru(bpy)3^*.  As  noted 
above,  the  radiative  and  nonradiative  processes  affecting  the 
lifetime  of  Ru(bpy)3^*  are  dependent  on  the  medium  in  which 
Ru  ( bpy )  3^^  res  ides . 

The  structural  information  gained  in  these  studies  can  be 
related  to  the  proposed  models  for  the  morphology  of  Nafion.  For 
example,  Gierke®^  developed  a  two-phase  model  (the  cluster- 
network  model)  of  the  morphology  of  Nafion;  according  to  this 
mode,  the  ionic  clusters  are  spherical  (ca.  40  A  in  diameter)  and 
are  interconnected  by  narrow  channels  (ca.  10  A  in  diameter) . 

The  cluster-network  model  assumes  that  both  the  ions  and  the 
solvent  are  localized  within  the  clusters  and  channels. 

Yeager  and  Steck®’  proposed  a  three-phase  model  where  a 
region  of  fluorocarbon  backbone  material  (with  some 
microcrystallites)  coexists  with  an  ionic  cluster  region 
containing  the  majority  of  the  -SO3'  sites,  counter  ions  and 
sorbed  water.  An  interface  is  presumed  to  exist  between  these 
two  regions.  This  interface  is  an  amorphous  fluorocarbon  region 
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containing  pendant  side  chain  material,  some  water,  and  some 
sulfonate  sites  with  associated  counterions.  The  clusters  are  no 
longer  assumed  to  be  spherical  and  intrusion  of  the  fluorocarbon 
phase  into  the  ionic  clusters  occurs.  Our  data  support  this 
interfacial  model  where  there  is  a  gradual  transition  between 
purely  fluorocarbon  and  purely  aqueous  environments,  it  is  this 
gradual  transition  which  produces  the  continuous  distribution  of 
chemical  microenvironments  within  the  polymer  phase. 

These  findings  have  some  relevance  to  the  electrochemical 
response  of  Ru(bpy)3^*  in  Nafion.  The  results  of  this  paper 
suggest  that  Ru(bpy)3^*  resides  in  a  variety  of  distinct  chemical 
microenvironments  within  Nafion  membrane.  The  energetics  of  the 
redox  reaction  are  clearly  subject  to  the  influence  of  these 
distinct  microenvironments  present  in  Nafion.  For  example,  the 
redox  waves  can  become  broaded  by  the  distribution  of  activation 
free  energies  for  the  redox  reaction.®®  While  broadening  of  the 
Ru(bpy)3^^  wave  in  Nafion  is,  indeed,  observed,  the  chemical 
microenvironments  are  not  so  distinct  as  to  produce  distinct 
redox  waves  or  shoulders. 
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TABLE  I.  Kinetic  and  Statistical  Parcuneters  Obtained  by  Fitting 
a  Biexponential  Decay  Model  to  the  Experimental  Ru(bpy)3^^ 
Emission  Data. 


a 

Ru(bpy)3^*,% 

^  short 

(%)  / 

b 

ns 

^  long  i^)  r 

b 

ns 

c 

2 

X 

0.01 

42.8 

+ 

11 

(2.75) 

442 

+ 

20 

(97.25) 

0.993 

0.03 

43.7 

+ 

13 

(1.81) 

475 

+ 

9 

(98.19) 

1.237 

0.05 

127 

+ 

48 

(6.07) 

505 

+ 

12 

(93.93) 

1.003 

0.07 

219 

+ 

53 

(16.20) 

613 

+ 

27 

(83.80) 

1.266 

0.10 

222 

+ 

42 

(17.87) 

633 

+ 

25 

(82.13) 

1.389 

0.20 

201 

+ 

24 

(15.25) 

606 

+ 

15 

(84.75) 

1.538 

0.30 

217 

+ 

27 

(16.73) 

608 

+ 

15 

(83.27) 

1.724 

0.50 

284 

+ 

26 

(32.34) 

758 

+ 

36 

(67.66) 

2.051 

0.80 

259 

+ 

24 

(23.85) 

657 

+ 

18 

(76.15) 

2.174 

®  Percent  of  -SOj' sites  in  Nafion  occupied  by  Ru(bpy)3^‘^. 

^  Error  in  lifetime  value  given  as  one  standard  deviation.  Value 
in  parenthesis  is  the  percent  of  that  lifetime  component  to  the 
total  luminescence  decay. 


Reduced  chi  square  value. 


TABLE  II.  Dependence  of  Emission  on  the  Ru(bpy)3^^  Loading 
Level  in  Nafion 


a 

Ru(bpy)3^^,  % 

•^max  r 

0.01 

593 

0.03 

597 

0.05 

599 

0.07 

600 

0.10 

600 

0.20 

602 

0.30 

601 

0.50 

602 

0.80 

603 

Percent  of  -SOs' sites  in  Nafion  occupied  by  Ru(bpy)3^*. 


TABLE  III.  Kinetic  and  Statistical  Parameters  Obtained  by 
Fitting  the  Dispersed  Kinetics  Model  to  the  Experimental 
Ru(bpy)3^^  Emission  Data. 


Ru(bpy)3^*, 

a 

% 

It  X  10®,  s 

-1 

T  / 

ns 

1 

b 

2 

X 

0.01 

2.639 

379 

+ 

18 

1.226 

1.003 

0.03 

2.241 

446 

+ 

8 

0.842 

1.204 

0.05 

2.365 

423 

+ 

6 

0.882 

1.002 

0.07 

2.264 

442 

+ 

5 

0.826 

1.262 

0.10 

2.292 

436 

+ 

4 

0.862 

1.384 

0.20 

2.329 

429 

+ 

3 

0.894 

1.538 

0.30 

2.306 

434 

+ 

2 

0.826 

1.718 

0.50 

2.332 

429 

+ 

2 

0.740 

2.066 

0.80 

2.278 

439 

+ 

2 

0.735 

2.189 

“  Percent  of  -SO3'  sites  in  Nafion  occupied  by  Ru(bpy)3^^. 
Reduced  chi  square  value. 


b 


TABLE  IV.  Kinetic  and  Statistical  Parameters  Obtained  by  Fitting 


the  Dispersed  Kinetics  Model  to  the  Experimental  RuCbpy),^^ 
Emission  Data  in  Oxygen-Free  Nation  Membranes. 


a 

Ru(bpy)3^*,  %  R 

X  10®, 

S'^  r  , 

ns 

7 

b 

2 

X 

0.01 

1.726 

579 

+ 

16 

2.157 

1.919 

0.03 

1.485 

673 

+ 

6 

0.919 

3.208 

0.05 

1.361 

735 

+ 

6 

0.764 

4 . 040 

0.07 

1.368 

731 

+ 

6 

0.868 

4 .208 

0.10 

1.372 

729 

+ 

11 

0.700 

2 .202 

0.20 

1.409 

710 

+ 

7 

0.619 

2.728 

0.30 

1.377 

726 

+ 

8 

0.675 

2.876 

0.50 

1.346 

743 

+ 

7 

0.633 

3.371 

0.80 

1.324 

755 

± 

7 

0.502 

3.707 

®  Percent  of  -SO3 

'  sites 

in  Nafion 

occupied 

by  Ru(bpy)3^\ 

b 


Reduced  chi  square  value 


Figure  1 .  Points-  Experimental  luminescence  decay  data  for 
Ru(bpy)3^*  from  Nafion  membrane.  Loading  levels  were  (A)  1  - 

0.80%,  2  -  0.50%,  3  -  0.30%,  4  -  0.20%;  (B)  1  -  0.10%,  2  - 
0.07%,  3  -  0.05%,  4  -  0.03%.  The  solid  curves  are  best  fit 
curves  obtained  from  the  biexponential  decay  model . 

Figure  2 .  Points  -  Experimental  luminescence  decay  data  for 
Ru(bpy)3^^  from  Nafion  membrane.  Loading  levels  are  (A)  0.80%, 
(B)  0.10%,  (C)  0.03%.  The  solid  curves  are  best  fit  curves 
obtained  from  the  Albery's  dispersed  kinetics  model. 

Figure  3 .  (A)  Luminescence  decay  curves  for  a  deaerated 

Ru(bpy)3^*  exchanged  Nafion  membrane.  (B)  luminescence  decay 
curve  for  an  aerated  Nafion  membrane. 
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APPENDIX 


DERIVATION  OF  ALBERY'S  EQUATION 

The  Arrhenius  expression  for  the  rate  constant  is 

k  =  Ae'^*^®^  (A.l) 

where  E,  is  the  activation  energy  and  A  is  the  preexponential 
factor.  Solving  for  Eg  gives 

Eg  =  -RTlnk  +  RTlnA  (A. 2) 

Assume  a  distribution  of  the  energies  of  activation,  the 
difference  between  an  individual  energy  of  activation  (Egi)  and 
the  average  energy  of  activation  (Eg)  is; 

Egi  -  Eg  =  -RTlnkt  +  RTlnA  -  (-RTlnk  +  RTlnA)  (A.  3) 

(E,i  -  Eg)/RT  =  -(Ink  -  Ink)  (A. 4) 

Now,  assume  that  the  energies  of  activation  are  distributed 
in  a  Gaussian  fashion.  Each  individual  energy  of  activation  can 
be  expressed  as  a  function  of  the  average  energy  of  activation 
and  the  standard  deviation  of  the  distribution,  7 


(Egi  -  Eg)/RT  =  -(7X) 


(A. 5) 


where  x  is  a  scaling  factor  (-«>  <  x  <  +<*>).  From  Equations  A. 4 
and  A. 5  we  obtain 

(7X)  =  (Inki  -  Ink) 
ki  =  ke^^ 

Equation  A. 7  will  be  used  later  below. 

Equation  A. 5  relates  each  individual  energy  of  activation  to 
the  average  energy  of  activation  and  7.  We  can  rearrange 
equation  A. 5  to  give 


(A. 6) 
(A. 7) 


X  =  (E.  -  E,i)/7RT  (A.  8) 

Now  let  us  assume  that  at  time  zero  we  have  this  Gaussian 
distribution  of  species  where  each  member  of  the  distribution 
decays  with  first  order  kinetics.  The  total  concentration  of 
these  species  at  time  zero  (0^(0))  is  given  by: 

C,(0)  =  C(0)exp(-x^)  (A. 9) 

Equation  A. 9  describes  a  Gaussian  where  C(0)  is  the  middle  of  the 
distribution  and  x  the  parameter  that  describes  the  distribution. 
The  first  order  rate  constant  for  all  the  species  with  a 
particular  x  value  is  given  by  equation  A. 7,  where  k  is  the  value 
at  the  middle  of  the  distribution.  To  obtain  the  equation 
describing  the  decay  for  the  whole  population  of  species,  having 
all  possible  x  values,  we  must  start  with  just  one  value  of  x. 


and  describe  the  first  order  kinetics. 

The  decay  law  for  species  corresponding  to  a  single  value  of 
X  follows  the  development  of  the  standard  first-order  decay  law 

dC^dt  =  -k^C,,  (A.  10) 


substituting  the  value  of  k^  from  equation  A. 7  we  obtain: 

dC^dt  =  -C,lce^  (A.  11) 

We  define  now  a  dimensionless  time  r  =  let  and  then  dt  =  dr/k. 
Equation  A. 11  becomes 


dC^(dr/k)  =  -C,lce^ 

Therefore 

dC^C^  =  -e^  dr 


(A. 12) 

(A. 13) 


Integrating  the  Ihs  of  equation  A.  13  from  C, (0)  to  C, (r)  and  the 
rhs  from  r  =  0  to  r  to  we  obtain 


and 


ln(C,(r)/C,(0))  =  -reT^ 
Cx(r)  =  Cx(0)  exp(-Te"'^) 


(A. 14) 

(A. 15) 


which  gives  the  concentration  of  species,  having  one  specific 
value  of  X  (i.e.  one  value  of  E^) ,  left  at  time  r.  To  obtain  the 
total  concentration  of  species  left  at  time  r  we  must  integrate 


4 


over  all  possible  values  of  x 


/fjKiO 

C(r)  =  r  Cjj(t)  dx=  I  CjcCO)  exp(-re'’^)  dx  (A.  16) 

J  -CO  J  -00 

We  know  that  0^(0)  =  C(0)exp(-x^)  (equation  A. 9)  and  therefore 

p-too  _  p-KO 

C(r)=  Cjj(r)dx  =  C(0)  exp(-x^)  exp(-re”’^)  dx  (A.  17) 

J  -00  J  -® 

C(r)  is  the  total  concentration  of  species  left  after  time  r. 
Rearranging 


C(r)/C(0)  =  J  exp(-x^)  exp(-re^)  dx 


(A. 18) 


The  total  concentration  of  species  at  time  =  0  is 


C(0)  =  J"*  CjjCO)  dx 


(A. 19) 


and  substituting  the  distribution  function  of  equation  A. 9  in 
equation  A. 19  we  obtain 

C(0)  ==  C(0)  J  exp(-x^)  dx  (A. 20) 

Rearranging 


C(0) 


C(0) 


exp(--x^)  dx 


(A. 21) 


Using  equation  A. 21  we  have  that 


C(r)/C(0) 


[C(r)/C(0) ] 


exp(-x^)  dx 


(A. 22) 


and 


C(r)/C(0) 


C(r)/C(0) 


exp(-x^)  dx 


(A. 23) 


Substituting  equation  A. 18  in  the  numerator  of  equation  A. 23  we 
obtain  Albery's  equation: 

exp(-x  )  exp[-r  exp(7X)]  dx 

C/C,  =  - - - -  (A.  24) 

f  exp(-x^)  dx 


where 


f  exp(-x^)  dx  = 
J  •CO 


Note  that  by  assuming  a  Gaussian  distribution  in  In  k  only 
two  adjustable  parameters  are  needed  to  describe  the 
distribution:  k,  the  rate  constant  of  maximum  probability,  and  y, 
the  measure  of  the  width  of  this  distribution.  As  noted  in  the 
experimental  section,  these  parameters  were  obtained  by  fitting 
simulated  and  experimental  transients. 


